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BERNARD SHIFFMAN, TATSUYA TATE, AND STEVE ZELDITCH 

Abstract. We determine the asymptotics of the joint eigenfunctions of the torus action on 
a toric Kahler variety. Such varieties are models of completely integrable systems in complex 
geometry. We first determine the pointwise asymptotics of the eigenfunctions, which show 
that they behave like Gaussians centered at the corresponding classical torus. We then show 
that there is a universal Gaussian scaling limit of the distribution function near its center. 
We also determine the limit distribution for the tails of the eigenfunctions on large length 
scales. These are not universal but depend on the global geometry of the toric variety and 
in particular on the details of the exponential decay of the eigenfunctions away from the 
classically allowed set. 
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1. Introduction 

A problem of considerable interest in both mathematics and physics is to determine the 
asymptotics of the distribution functions 

Dj(t) := Vo\{z : \<Pj(z)\ 2 > t} 
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of an orthonormal basis of eigenfunctions of a Laplacian (or similar Hamiltonian) on a 
compact manifold (see pp). In general, it is hopelessly difficult to obtain more than crude 
bounds on such distribution functions, which of course control the £ p -norms of the eigen- 
functions. Numerical analyses and heuristics from quantum chaos and disordered systems 
suggest however a rich picture in which the asymptotics of Dj(t) is related to the classical 
dynamics underlying the eigenvalue problem. (Some references will be discussed at the end 
of the introduction.) The purpose of this paper is to give a rather complete analysis of the 
limit distribution of eigenfunctions in one of the few settings where such a detailed analysis 
is possible, namely where the phase space is a toric Kahler variety (M,u). 

Let us recall the definitions (see §1 for details). Toric varieties are complex manifolds on 
which the complex torus (C*) m acts with an open dense orbit. By a toric Kahler variety 
we mean a toric variety equipped with a Kahler form (M, uj) that is invariant under the 
underlying real torus T m . The action of T m is Hamiltonian with respect to u, and thus toric 
varieties are models of completely integrable systems. They are of a very special type because 
integrable systems usually generate an IR m action rather than a T m action. Although there 
are rigidity theorems limiting the class of such examples in the world of real Riemannian 
manifolds |LSj . toric varieties provide a plentiful collection in the world of complex manifolds. 

The torus action can be 'quantized' or linearized on the Hilbert space completion of the 
coordinate ring 

oo 

H:=($H°(M,L N ), (1) 

N=0 

where L — > M is a holomorphic line bundle with C\(L) = and where H°(M, L N ) denotes 
the space of holomorphic sections of its iV-th tensor power. This quantization is generally 
known as the holomorphic (Bargmann-Fock) representation in the physics literature. The 
space 7i is spanned by joint eigenfunctions of the linearized (C*) m action, which we refer 
to as 'monomials.' In the fundamental case of M — CP m , the joint eigenfunctions are the 
monomials given in an affine chart by 

Xa :C m ^C, Xa{z) = z a . (2) 

The monomials lift (by homogenization) to homogeneous monomials on C m+1 . 

We consider the case where M is a smooth projective toric variety; i.e., M = Mp, L = Lp, 
where P is an integral Delzant polytope (see Then the linearized T m action is generated 
by m commuting operators Ij, j = 1, ... ,m on Mp which preserve H°(Mp, Lp), and the 
joint spectrum of the eigenvalue problem 

^« = <¥P«> « = (ai,... a m )eR m , <f^eH°(Mp,L^) (3) 

consists of lattice points a G NP fl Z m . 

Our main results concern the asymptotics of their distribution functions 

£> 7 (t) :=Vo\{zeM P : |<^)| 2 > t} (4) 

with 7 G NP fl Z m as iV — > oo. The function |^(2;)| 2 is often called the 'Husimi distri- 
bution' in the physics literature, and thus our results determine its distribution law. The 
norm ^^(z)! of (p~(z) G Lp is the pull-back of the Fubini-Study norm under a monomial 
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embedding of the form 

$P = [Ca(l)Xa(l), ■ ■ ■ , C a{d+1)X a( d+ l)} : (C*) m -> CP d , P n TT = {a(l), . . . , a(d + 1)} , 

(5) 

for a choice of constants c a (j) G C*. The volume in Mp and the Hermitian norm hp on Lp 
are by definition the pull-backs of the Fubini-Study metric and form under this monomial 
embedding, and the C 2 norm on the space H°(M,Lp) is in turn induced from the volume 
form and the Hermitian pointwise norm of hp. (See £|21for details.) 

As Figure ^ illustrates, the monomial ^ is something like a Gaussian bump centered on 
the real torus fXp 1 (a), where /ip : Mp — > P is the moment map for the classical Hamiltonian 
T m -action on M P (see §1). 




Figure 1. 2-k/N times the monomial |^ Q (^)| 2 for iV = 1 (left) and N = oo 
(right) for m = 2 and a = (2,3), where S is the standard simplex. The 
variable z is chosen as u t— > z = e (p<±+u/\/N)/2^ g ee p r0 p OS jtion 13. 15l in Section 
□ 



Since we would like to determine the properties of eigenfunctions (p a G NP when a is 
large, but not necessarily a multiple of a lattice point in P, we shall consider a sequence of 
approximate multiples, as in the following definition: 

Definition 1.1. Let on G NP fi Z m be a sequence of lattice points, and let x G P. We say 
that {oat} is a sequence of approximate multiples of x if 

a N = Nx + 0(l). (6) 

Our first result gives the pointwise behavior of the eigenfunctions: 

Theorem 1.2. Let x be a point in the interior P° of the polytope P, which is not necessarily a 
lattice point. Then there exists a non-negative function G C°°((C*) m ) such that b^(z) = 
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if and only if z G /ip 1 (x), and for every sequence E NP D Z m 0/ approximate multiples 
of x, we have 

1*4 (*)| 2 = c(P,x) (^J e- Nb ^ [1 + O(N-i)] 
uniformly on (C*) m , where c(P,x) G M + . 

The constant c(P, x) is denned in fHHj) and (jlo]) (and also appears in Theorems 11.31 and 
OD- 

When we consider a sequence of the lattice points of the form = Na with a lattice 
point a, the assumption that a G P° is not necessary. In fact, we will give the pointwise 
asymptotics for a^ = Na with any lattice point a G P fl Z m in Section (see Propositions 
E2andEEJ). 

In the above Theorem, the function b x (z) (x G P°) on (C*) m is defined by 

where t x (z) G M m is the vector given by the equation 

lj P (e- T * {z)/2 z) = x, 2 G (C*) m • (8) 

Here \ip is the moment map for the T m action on Mp (see (|23|) for the definition), and we 
write 

e r z = (e ri Z!, . . . , e rm z m ) for r G R m , z G (C*) m . (9) 

We can express (as in jBZ2, (17)— (18)]) the function b x in the more intuitive form as follows: 
we introduce the real power 'monomials' 

\Xx{z)\ ■= \z\ x = \z x \ xi ■■■\z m \ Xm 

and define 

M P x {z) := Wx } Z) ) P M , (10) 
sup \ Xx {z)\ P 

\ Xx {z)\ P := \X^}L= (ze(C*) m ). 

(The normalized monomial A4 P has sup-norm 1, attained on the torus fj,p (x).) Then ((7j) is 
equivalent to: 

b£(z) = -2]ogM* M (z), (11) 
Since the sequence of monomials flattens out exponentially quickly away from the peak 
set /ip 1 (x), the distribution function is clearly tending to zero. The rate of decay of the 
distribution function is given by the following result: 

Theorem 1.3. (i) Let a^ G NPdl/ 71 be a sequence of lattice points which are approx- 
imate multiples of x G P° (see ©j. Then, fort > 0, we have 

(vrm)" 1 / 2 (\ogN\ m/2 

Da N {t) 



where (cf. (|2fl|l ) 



c(P,x)T(m/2 + l) V N 
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(ii) Let aedPnZ m . Then, for t > 0, we have 

(nd(a)) d ^/ 2 



D Na (t) 



logiV 



d(a)/2 



c(P,a)T(d(a)/2 + 1) V N 

where d(a) := m + codimF a , F a being the face of P containing a. 

Here ~ means the ratio of the left and right hand sides tends to 1, and the constants 
c(P, x), c(P, a) are given by (J4"3~j) and fffOj) — . We recall that if x is a point in a face F of 
codimension r, then /ip 1 (a;) = T m_r . Hence d(a) = 2m — dim/ip^a). 

The exponentially localized behavior of the monomials suggests studying the distribution 
function on various length scales. First, we show that the D aN have a universal scaling limit 
on a small length scale: 

Theorem 1.4. (i) Let cin G NP n Z m be a sequence of lattice points satisfying the 
condition © with some point x G P°. Then, for < t < c(P,x), we have 

lim {N/27r) m/2 D aN ({N/27T) m/2 t) = — — 

N— >oo 

(ii) Let a G Pf]Z m . Then 



c(P,x)T(m/2 + l) 



{\og{c{P,x)/t)T/\ 



lim (N/2tc 



.d(a)/2 



D 



Na 



(N/271 



d(a)/2 



1 



(log(c(P,a)/t)f a)/2 



n^oc" ' ' — \^ ' ' / c(P,a)r(d(a)/2 + l) 

for < t < c(P,a), where d(a) = codim/ip (a). 
A sample graph of the scaling limit distribution function for P = 7S is given in Figure 121 




Figure 2. Scaling limit distribution for |y?^J 2 with m = 2, a = (2, 3), P = 7S 

As in Theorem II .41 the limit of the rescaled distributions has a universal form, i.e. it does 
not depend on the geometry of the manifold Mp, and is given by a logarithmic power of 
the form (logc/t)^ 2 with some constant d. The logarithmic power appears because the C 2 - 
normalized monomials are close to Gaussian around the peak set /ip X (x) on a vector space of 
dimension m/2 (or dimension d(a) for the boundary lattice cases). More precisely, Theorem 
11.21 (and Propositions 13. 51 and I3.8JI shows that the function has a positive definite Hessian 
at the peak point. We then observe that the distribution of a Gaussian function on M. d , 

p -{Au,u)/2 
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(where A is a real positive (d x d)-matrix) is given by the logarithmic power law 



v A (u G M d ; g(u) > t) = — — flog -) 

v ' yv ' ' cT(d/2 + l) V B tJ 



c = , < t < c, 

VdetA 

relative to the normalized Lebesgue measure 

det A , 



VA : (2vr) rf / 2 



Thus, the rescaled distribution of an £ 2 -normalized monomial at the 'center' of its localized 
bump has a universal Gaussian form. 

To analyze the 'tails' of the eigenf unctions, we next use an exponential rescaling of the 
distribution function so that the global distribution law has a non-zero limit as iV — > oo. As 
may be expected, it is no longer universal but depends on the geometry of (M P ,u P ). 

Theorem 1.5. Let a N E NP n Z m satisfy the condition (0) with a point x G P°. Then 



where A{p) is the Hessian matrix oflog^^^p \cp\ 2 e^ ,p ' . 

A more general scaling limit law is given in Theorem 14.31 In the above theorem, the assump- 
tion that x G P° is not necessary. In fact, we will give similar result for = Na with any 
lattice point a 6 P ( Theorem 14. 3j) in Section 0] 

Our strategy for proving Theorems II . 31 and II . 51 on the distribution functions of monomials 
is based on their pointwise asymptotics (Theorem 11.21 and also Propositions 13.51 and 13.81 in 
Section El). Pointwise asymptotics of monomials are more or less equivalent to asymptotics 
of their L 2k norms. Since the latter are of independent interest, we state these asymptotics 
explicitly in: 

Theorem 1.6. (i) Let a N E NP n Z m satisfy the condition (jUJ) for a point x E P°. 
Let \\(p^\\2k denote the C 2h -norm of the C 2 -normalized monomial (p^ with the weight 
7 G NP. Then we have 



We close the introduction with some general remarks and references. As mentioned above, 
our results pertain to phase space distribution of eigenfunctions (Husimi distributions) rather 
than to their configuration space distribution. To our knowledge, the only prior example 
in dimension > 1 for which the limit distribution of eigenfunctions has been determined is 
the case of certain (so-called) Hecke eigenfunctions of discrete quantum cat maps, due to 




det A(p) dp , 




where Ok{N 1 ) depends on k. 
(ii) Let a G P fl Z m with d(a) = codim/Xp 1 (a). Then we have 
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Kurlberg-Rudnick |KRj . They work in a simpler discrete model rather than the holomor- 
phic model. The main result of Kurlberg-Rudnick |KRj is that the distribution functions 
of (un-scaled) eigenfunctions tend to the semi-circle law. Their method was to relate the 
eigenfunctions to exponential sums studied by Katz |Kaj and to apply the value distribu- 
tion of exponential sums. Since value distribution depends on the representation, it is not 
clear that the same semi-circle law would hold for Hecke eigenfunctions in the holomorphic 
(Bargmann-Fock) representation, and this appears to be a challenging and interesting prob- 
lem. Part of the motivation for this paper was to set a baseline for eigenfunction distribution 
problems by studying a class of explicitly solvable examples. 

It would also be interesting to study the limit distribution of real eigenfunctions in the 
Schrodinger representation, i.e. on the configuration space rather than the phase space. 
To our knowledge, the physics results mainly pertain to these configuration space results. 
The cases most studied and speculated about are those of chaotic or disordered systems. 
When eigenfunctions are delocalized, their spatial distribution is conjectured to be Gaussian 
(see e.g. |Be| IFE| IHel IHR| IMi| ISS[ IPAj ) . In the opposite regime where the eigenfunctions of 
disordered systems are exponentially localized, the expected distribution in the low amplitude 
(tail) region is given by a power of a logarithm JVIFJ, precisely the one we obtained in the 
high amplitude (center) region. The reason is that the distribution law is universal and 
Gaussian in the tail region for exponentially localized eigenfunctions of disordered systems 
|MFj . while it is universal and Gaussian in the center for our problem ( Theorem 1 1.2|) . 

The only studies we have located which are related to distribution laws of integrable 
eigenfunctions are those of one of the authors with J. A. Toth (cf. |TZj ) and that of Berry- 
Hannay-Ozorio de Almeida (BHQ j , which describes the asymptotic expansions of C 2p norms 
(moment intensities) of real oscillatory integrals of several stable types. Such oscillatory 
integrals define quasimodes for a quantum integrable system, and in generic cases one can 
express eigenfunctions as oscillatory integrals of various kinds 1TZ' . In general, many possible 
kinds of oscillatory integrals could arise, including ones associated to singular Lagrangean 
tori. Hence, one can only expect complete results in special cases. To take the simplest 
example, our methods could be adapted to find the scaling limit distributions of squares of 
the standard spherical harmonics on S 2 (or S m ). These eigenfunctions are also joint 
eigenfunctions of commuting operators which generate a quantum torus action. To our 
knowledge, the distribution laws of even such simple eigenfunctions are unknown at this 
time. 

Acknowledgments: This paper was written during a stay of the second author at Johns 
Hopkins University on a JSPS fellowship. He would like to express his special thanks to the 
faculty in the Department of Mathematics of Johns Hopkins University. 

2. Background on toric varieties and moment polytopes 

We summarize here some basic facts and terminology on toric varieties from |STZj . Recall 
that a toric variety is a complex algebraic variety M containing the complex torus 

(CT := (C \ {0}) x • • • x (C \ {()}) 

as a Zariski-dense open set such that the group action of (C*) m on itself extends to M. 
We consider here smooth projective toric varieties; they can be given the structure of a 
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symplectic manifold such that the restriction of the action to the underlying real torus 

T m = {(0, . . . Cm) G (CT : 101 = 1. 1 < 3 < W 
is a Hamiltonian action (see §2.2)1 . These toric varieties can be constructed from Delzant 
polytopes either by symplectic reduction (see |Gu| ) or by gluing affine toric varieties described 
by the normal fan of the polytope (see jFuj ) . However, for our analysis, it is more convenient 
to define the toric variety Mp associated to the Delzant polytope P through a monomial 
embedding as follows (see GKZ , Chapter 5]). Suppose that P is a Delzant polytope and let 

P n Z m = {a(l), a(2), . . . a(d + 1)} . 

We shall write J2 a eP = Saepnz™- (Recall that a Delzant polytope is a convex integral 
polytope in M. m with the property that each vertex is incident to exactly m edges and the 
primitive vectors in Z m parallel to these edges generate Z m .) 

To define the monomial embedding, we fix an arbitrary c = (ca^^^a+i)) G (C*) d+1 . 
Then, we define 

$ P = [Ca(l)Xa(l), • • • , Ca(*fl)X«(*fl)] = (C*) W - CP d | (12) 

i.e., 

d> P (z) = [c a{1) z a ^,...,c a{d+1) z^ d+ \ z e (C*) m - 
The toric variety Mp = Mp is defined as the Zariski-closure of the image $p((C*) m ) of the 
monomial embedding $p in the complex projective space CP d . 

Since our polytope P is assumed to be Delzant, $p is an embedding and the variety Mp 
is smooth. The symplectic (or Kahler) form on Mp is given by 

up = up := $pW FS , (13) 

where c^fs = 2~9d log ||C|| 2 denotes the Fubini-Study Kahler form on CP d with homogeneous 
coordinates (Co, ■ ■ ■ , Cd)- On (C*)' m , we have 



^p = ^^log^|c a | 2 |^| 2 . (14) 

aeP 

The volume form on Mp is given by gTVoIm = ^w™. 

2.1. The line bundle Lp and associated circle bundle Xp. We define the line bundle 
Lp — > Mp by Lp = Lp := $p*C(l), where C(l) denotes the hyperplane section bundle 
on CP d . Recall that the space of holomorphic sections H °(CP d , 0(1)) consists of the linear 
functions A : C d+1 — > C, and that the Fubini-study metric on 0(1) is given by 



|a!fs([CD = W (C g c d+1 ; 



which has curvature form cups = T^:<9<91og ||CI| 2 - We endow Lp with the Hermitian metric 
hp := $p/iFS, which has curvature up. 

Each monomial Xa with a e P D Z m corresponds to a section of H°(Mp, Lp) and vice 
versa. To explicitly define this correspondence, we make the identifications: 

X P a{j) = c-^Q G H°(M P , Up) = &*H°(CF d , 0(1)) , 1 < j < d + 1 . (15) 

More generally, 

$p*#°(CP d , 0(iV)) ^ H°(M P ,L^) (N>1), (16) 
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and a basis for H°(Mp, Lp) is given by the sections {x^ '■ 7 G NP R Z m } corresponding to 
the monomials {x-y}- These sections are given by 

where ft,...,^GPnZ m such that 7 = + ■ ■ ■ + (3 N (seejEH RT7] ). 

So far, we have not specified the constants c a . For studying our phenomena, the choice of 
constants defining the toric variety Mp is not important. However, when our polytope P is 
the full simplex pS, we shall use the special choice c a = (^) 1//2 , where ( p ) is the multinomial 
coefficient (see H2A\f . 

The associated principal 5 1 -bundle Xp = Xp of the line bundle Lp — > Mp is defined by 

^^{(^JeLp 1 ; Hp = 1}, 

where |u| P denotes the norm of t> with respect to the Hermitian metric on Lp 1 induced by 
hp. 

We identify sections sn of with equivariant functions sjv on X by the rule 

Sjv(a) = (A®^, s N (zj) , \eX z . (17) 

Clearly, s N (e ie • x) = e iNe s N (x) if s N G H°(Mp,L^). It should be noted that for each 
s G H°(Mp, Lp), we have 

|sjvOe)| = l s iv(^)|p 

where a; G Xp is in the fiber over z G Mp, |sat(^)|p denotes the norm with respect to the 
Hermitian metric on Lp induced by the metric hp. 

In particular, for a G P (7 Z m , the 'monomial' Xa £ H°(Mp, Lp) given by lifts to an 
equivariant function on the circle bundle Xp — > Mp, and we write 

:= C ^U)Xa{j) = ^ (18) 

where t P : Xp -> is the lift of the embedding Mp <^-> CP ' (d = #P - 1), that is, t P is 

the restriction to Xp of the natural inclusion Lj, 1 ^ 0(— 1). (Of course, depends on c, 
which we omit to simplify notation.) We also consider the monomials 

p p 

so that m p is the equivariant lift of m a to Xp. In terms of local coordinates (z, 9) on 
Tj- 1 {{C*) m ) C Xp, we have 

= - ° - 1/2 ■ (19) 

E^pI^I 2 ) 

Since its absolute value is independent of 0, we can write 



<{z)\ P = \m p a {z)\ = lCaZ 1 . (20) 

012 



1/2 • 



We give H°(M P , Lp) the inner product 

(si,s 2 >= / U(z),~^(z)) dVo\ M (z), Sl ,s 2 eH\M P ,L N p) , (21) 
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and the C? norm ||s|| = (s, s). We note that the sections 

{ X P a E H°(M P , Lp) : a = (a 1( ...,a m ) e NP} 

are orthogonal but not normalized. We normalize them to obtain an orthonormal basis for 
H°(M P , Lp) consisting of the sections 

rf = & ' (22) 

Their equivariant lifts (p^ form an orthonormal set of monomials on Xp. 

We include below a table of notation to help the reader keep track of the various monomials: 





monomials on C m 


sections of Lp 


monomials on Xp 


ae NP 


Xa(z) = Z a 


Xa 


Xa 


aeP 




m a = C »Xa ( N = 1) 


-~ P 
m a 


ae NP 




= Xa/\\Xa\\ 


ft 



2.2. Moment maps and torus actions. The group (C*) m acts on Mp and the subgroup 
T m acts in a Hamiltonian fashion. Let us recall the formula for its moment map fip : Mp — > 
]R m : on the open orbit (C*) m , we have 

Pp(z) = n%(z) = = ^-— 2 K\ 2 \z a \ 2 a = K(z)\ 2 a . (23) 

For any c, the image of Mp under up equals P. The moment map //p is invariant under the 
T m -action on Mp. By the identification (C*) m = T m x W m , the moment map fip defines a 
map from M m to P, which is a diffeomorphism between M m and the interior P° of P. 

The action of the real torus T m lifts from Mp to Xp and combines with the S 1 action to 
define a T m+1 action on Xp. Recall that under the monomial embedding 

$ P : (CT ^ M P - CP d , * i ► [c o(1) * a W . . . , c a(d+1) ^ d+1 )] , 

the T m action on Mp C CP d is given by 

■ [Ci, • • • , Cd+i] = [e i(a(1) '^Ci, • • • , e^^^Ofi] • (24) 
The action (|24|) lifts to an action on Lp 1 : 

■ C = (e^W'^Ci, • • • , e i<Q((m) ' v> Cm) • (25) 
Since the circle bundle Xp C S +1 is invariant under this action, (|25|) also gives a lift of the 
action J2D to Xp. 

We also have the standard circle action on Xp: 

e id ■( = e w ( , (26) 

which commutes with the T m -action (|25|). Combining (|25|) and (|26|). we then obtain a 
T m+1 -action on Xp: 

(e ie ,e i(p \...,e i,fim )»C = e ie (e^ -C) ■ (27) 
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2.3. Fourier analysis. In this section, we shall explain an aspect of Fourier analysis on 
toric varieties which describe the complete integrability of the system. 

The Hardy space 7i 2 (Xp) is the Hilbert space spanned by the equivariant lifts of sections 
to functions on Xp with the inner product 

(f,9) = [ f9dV, dV = a P A (da P ) n -\ 
J x 

where ap is a contact 1-form defined by the Hermitian connection on Lj, 1 such that dap = 
Ti*uj P . Under the identification TC = H 2 (Xp), the inner product is the same as the in- 
ner product on H = ^ N H°(M Pl L^) given by (ED)- Alternately, H 2 {Xp) consists of the 
functions F E C 2 (Xp) satisfying d b F = (see e.g., [SzTUZe]). 

Under the S 1 action, the Hardy space then has the orthogonal decomposition 

oo 

n 2 (X P ) = Q)H 2 N (Xp), (28) 

N=0 

where TC%(Xp) consists of elements s E TC 2 (Xp) such that s(e t9 -x) = e tN0 s(x). We recall that 
the (equivariant) 'Szego projectors' are the orthogonal projection onto H°(Mp, L c p) = 
n 2 N (Xp). If {5f } denotes an orthonormal basis of H°(Mp, Lf), and Sf denote their lifts 
to X, then the projector is given by the kernel function 

U N (x,y) = J2sf(z)Sf(y) : C 2 {X P ) - H 2 N (X P ). (29) 

3=1 

We now describe how one can combine the eigenvalue problems given by Q for varying N 
into a homogeneous scalar eigenvalue problem on X. To do this, we define for each N e N, 
the 'homogenization' NP C Z' m+1 of the lattice point in the polytope NP to be the set of 
all lattice point 5* of the form 

= a := (Np — |ot|, ai,... , a m ), a = (oti, . . . , a m ) E NP D Z m , 

where p > max^ g p n gm |/3|. We also define the cone Ap = U^=i NP. It is well known that 
rays Na in this cone define a semiclassical limit. 

In this section, we use the more precise notation (p&{x) for the £ 2 -normalized monomial 
<Pu(. x ) (since iV is not specified in the latter), for a E Ap. 

The torus action on Xp can be quantized to define an action of the torus as unitary 
operators on 7i 2 (Xp). Specifically, we let I±, . . . I m denote the differential operators on Xp 
generated by the T m action: 

(lMO = -^-S(e i<p - C)l^=o, SEC^(Xp). (30) 

We recall the following observation from [BTZ : 

Proposition 2.1. For 1 < j < m, 

(i) Ij : H 2 N (Xp) -> H 2 N (Xp); 

(ii) The lifted monomials tp? E TLjj(Xp) satisfy Ij(f>& = ctj<p& (o. E A P ). 
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Furthermore, we note that 

^:H 2 n (X c p )^H 2 n (Xp) , ]^0 §N = N§N for s N eU 2 N (X c P ) . (31) 
Thus, the monomials <p& are the solutions of the joint eigenvalue problem 

Ijfia = Ujfia, Ot e M m+1 , Bbfia = 0, j = 0, . . . , 171 (32) 

with the commuting operators: 

d m 

I° = -QQ~^2li> Il,...,I m - (33) 

The joint eigenvalues are the lattice points a G A P . 

2.4. Monomials on projective space. In the case of CP m , the £ 9 -norms of the monomials 
V^vi can be evaluated explicitly in an elementary way. We give the details in this section. 

When the polytope P is the unit simplex S, we have = CP m ; furthermore is the 
hyperplane bundle 0(1). We can identify L^ 1 = C r™{-l) with C m+1 with the origin blown 
up, and the circle bundle X s C L^ 1 is identified with the unit sphere S 2m+l C C m+1 . The 
equivariant lifts to X s of sections of 0(N) = consist of homogeneous polynomials 

nco, . . . , u = Yl c ^ x (c A = Co° • • • ct) 

\X\=N 

inm + 1 variables. The induced Fubini-Study metric on O(N) is given by 

|F(C)| s = |F(C)|/||Cir, for FeH (CF m ,O(N)). 

(Here the constants c a are taken to be 1.) Identifying F with the polynomial f(z) = 
F(l, zi, . . . , z m ), the norm can be written 

l/W| s = l/WI/(i + lkllT /2 (^q, 

and the inner product is given by: 



^ = i/ TTTW^^W, f,geH (CF m ,O(p)). (34) 

IW. J£m \L -f- \\Z\\ f 

We are interested here in the case P = pS. Then M p £ = CP m , and the line bundle L p s 
is identified with the p-th tensor power 0(p) of the hyperplane section bundle. The circle 
bundle X p s is the lens space X p s = S l2m+1 /Z p . By lifting equivariant functions from the 
lens space X pS = S 2m+1 /Z p to the sphere S 2m+1 , we see that H 2 N (X pll ) 2* H 2 Np {S 2m+1 ). 
Hence, we shall replace the iV-th Hardy space H 2 N (X P ^) by H 2 Np (S 2m+l ) below. Then the 
equivariant lift ■ S 2rn+1 -> C of x v * e #°(CP m , 0(p)) is given by the homogenization: 

X P a{x)=x a , a = (p - \a\,ax, . . . ,a m ). 

In this case, we shall use the special choice of the coefficients of the monomial embedding: 

p\ 2 ( p\ p\ 



\aj \a J (p — \a\)\ai\ ■ • ■ a m \ ' 
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so that by (JHJ), we have 



: ^io g (E, a ,< p Qkl 2 ) = ^r^Mi + \\4 2 f = ^fs- 

Furthermore, 

/i P sW :=< S (2) = ^ pr— p X, ( a ) |z " |2a= 1 + Vkp tN 2 ''^!^). (35) 

l^\a\<p\a)\ I | a |< p \ / ^' Jl 

where the last equality follows by differentiating the identity (1 + J2 x j) P = Yl\ a \<p (a)*^ - 
Note that this choice gives us the scaling formula 

As before, for each m-dimensional multi-index (3 with \/3\ < p, we define an (m + 1)- 
dimensional mult i- index (3 by 

An), A)=P-|/3|, & = & (j = l,...,m). 

Recall that |^ E (x)| = |^ S (^)| pE with ar G S 2m+1 , vr(x) = z G CP m . 

The £ 9 -norms of the monomial can be evaluated explicitly as follows. 

Proposition 2.2. For a e p£ n Z m , we /iai>e t/te precise formula: 

■(Np + m)\] q/2 nr =0 r(iVga,/2 + l) 



I 119 



(iVa) 



pm(q/2-i)T(Npq/2 + m + 1) ' 



Proof. We write £jvs = XjVa> which we consider as a function in 7-^ 2 Vp (S' 2m+1 ). For q > 1, we 
set 



!,(#)= / e-l^lxiva^)! 9 ^), 

JC m+1 

where <i£(x) denotes Lebesgue measure on C m+1 . We shall compute the integral T q {N) in 
two ways. First we use polar coordinates on C m+1 . The measure di is expressed as 

dt{x) = r 2m+1 dr da, x = ra, r > 0, o e S 2m+1 . 

Then we have 

1. 



g2m + l 



2 



r(iVpg/2 + m + l) / |£iVd(a)| 9 rfa 



tj2m+l 



To relate the volume form da on S 2m+1 to that on X P £, we recall that u p s = puips and 
therefore the volume form on X pE is p m times the Fubini-Study volume. Recalling our 
convention that Vol(X s ) = Vol(CP m ) = ^j, we then have dVo\x pTj = 27 fm+i < ^ cr ' anc ^ hence 

X *( N ) = ^t v ( n p<i/ 2 + m + VWxNaWl- 
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On the other hand, if we use polar coordinates for each component Xj of x G C m+1 , we 
have 

m „ 

T q (N) = TT / e~ lxl2 \x\ Najq d£(x) 
3=0 Jc 

m m 

= (27r) m+1 J| / e - r2 r Na i g+1 dr = vr m+1 JJ V{Na jq /2 + 1). 



Hence we obtain 



and therefore 



IC 119 — n m '■> 
\XN&\\ q —P 



UT =0 r(N qaj /2 + i) 



\XNc 



T(Npq/2 + m + l) 
2 p m (Na)\ 



(36) 



(iV> + m)V 

Since = XiVd/||XiVa||> the identity follows from (|3U j) - (|37j) . 

The following proposition is a direct consequence of Stirling's formula and Proposition 



(37) 
□ 



Proposition 2.3. Let a e n Z m , and let J = {j e Z : < j < m, atj ^ 0}, where 
(Xq — p — \a\. Set d(a) = 2m — #J. Then we have 



Wn 



a II q 



N \ 



- 



2tt J 



(q/2-l)d(a) p 9 /2-l/ 27r )» , (5-2) 

5/2=1' 



|^S||2 _ -J^P + H 
WNaWoo — P 



(Na)\ 

In particular, again by Stirling's formula, we have 



11 w 

pp 



N 



Wn 



p£ II 2 



(27T) 



P 



1/2 



2tt J 



d(a)/2 



In the next section, we will obtain similar formulas for monomials on general toric varieties. 



3. POINTWISE ASYMPTOTICS ON GENERAL TORIC VARIETIES 

We now consider the case of general toric varieties. Our first purpose is to find the 
pointwise asymptotics of the monomials and to prove Theorem II .21 We then asymptotically 
determine their £ 2fc -norms. 
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3.1. Pointwise asymptotics: interior points. First of all, we shall consider a sequence 
o/v of lattice points in NP. We assume that there exists a point x G P° such that 

a N /N = x + OiN- 1 ). (38) 

Under the assumption (|38p. the point a^/N is in the interior P° of the polytope P for 
sufficiently large N. Thus the analysis is performed on the open orbit (C*) m , and hence the 
coordinate 

z = e p/2+ilp e(C*) m , p, v eR rn 

will be useful. The moment map \ip is also invariant under the Hamiltonian T m -action, and 
it is well-known ( |Fuj ) that it induces a diffeomorphism: 

ftp : R m = (C*) m /T m -> P°, /ip(p) := ^ P {e p ' 2 ). (39) 
In these coordinates, the function 5^ defined by (J7J) can be written simply as 

b x (p) = f{x,p) ~ f(x,Px), Px=fip\x), f{x,p) : = logfc(p) - (p,x), (40) 
where the function is the 'polytope character' 

Kp)= E M 2e<p,/3> - ( 41 ) 

/3ePnZ m 

The vector t£(z) in (jHJ) is given by 

r^z) = p-p p x , z = (42) 
We also define the symmetric real m x m matrix 



2 



A(P,x) := \rfip(e^ /2 ) (3 ® (3 - x ® x. (43) 
/?eP 



Lemma 3.1. The real symmetric matrix 

A{p) := £ |mJ(e^ 2 )| 2 /3 <g> /3 - /x P (e^ 2 ) ® /x P (e p/2 ) 
/3eP 

is positive definite, for all z G (C*) m . 

Proof. We must show that (A (g> A, > for A G (M m )' \ {0}. Consider the vectors u, v G 

given by Mq = |m^(e^/ 2 )|, f Q = |m£(e"/ 2 )|A(a). Since ||w|| 2 = E ae p l™£ (e p/2 )| 2 = 1 
and /ip(e P//2 ) = XlaeP l^a (e P//2 )| 2 a, we have 

(A®A,A(p)) = ^|m:(e^)| 2 A(a) 2 -(^|m:(e^)| 2 A(a)) 

a£P VqGP / 

= |M| 2 -(m,v) 2 = ||u|| 2 ||<u|| 2 - («,<u> 2 > 0. (44) 

Since M a = |«^(e p / 2 )| 7^ for all a G P PI Z m and v a /u a = X(a) is not constant on P fl Z m , 
the Cauchy-Schwartz inequality in (jUj) is strict. □ 
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In particular, the matrix 

A{P,x) = A(p p x ) 

is positive definite. 

For a point x G P°, we now define the constant 

c{P,x) := . 1 (45) 
y/det A{P,x) 

Lemma 3.2. In the coordinates z = e p l 2+%e on (C*) m , the volume form on Mp is given by 

= T^w det A(p) dpdQ . 



Proof. We note that 

A(p) = '$2kf } (p)l3®l3- [J2h(p)A ® few/H = Hess p logA;(p), (46) 
/36P VfleP / \/3eP / 

where kp(p) = |m^(e P//2 )| 2 and k(p) = J^p^p |cs| 2 e( p,/3 ^. The conclusion follows from (|4b|). 
recalling that 

u P = $ C P ^ FS = ^dd log V |^i 2 | X/3 (^)i 2 = ^=5a91og*(p). (47) 
Z7T z — ' Zix 

□ 

It should be noted that b^(p) grows as \p\ — > oo, as stated in the following simple lemma. 

Lemma 3.3. Lei K d P° be a compact set. Then there exists positive constants R > 0, 
c> such that f(x, p) > c for (x, p) G K x M' m ; |p| > R. 

Proof. For any (x, p) G P° x M m , we define 

M(x,p) — max (p,3 — x). 

If x G P°, then the polytope P — x contains the origin in its interior. Thus, clearly we have 
M(x, p) > for any (x, p) G P° x (IR m \ 0). Next, we note that the function (x, p) h-> M(x, p) 
is continuous. To see this, let (x n , p n ) be a sequence such that (x n ,p n ) — » (a;, p) G P° x M m . 
Then, for any 3 G P n Z m , 

\(p n ,P-x n ) - (p,(3-x)\ < C\p n -p\ + \p\\x n -x\. 

By using this inequality, we can show that 

\M(x n ,p n ) - M(x,p)\ < C\p n - p| + \p\\Xn ~ x\. 

Now, for a compact set if C P°, we set 

M(K) = min M(z, p) > 0. 

(x,p)£KxK.™,\p\=l 



DISTRIBUTION LAWS FOR INTEGRABLE EIGENFUNCTIONS 17 

We set Co = min^gpnzm \c/3\ 2 . Since P fl Z m is a finite set, there exists (3 = /3(x,p) such that 
M(x,p) = (p,P(x,p) — x) for (x, p) with x G P° and |p| = 1. Thus, for x £ K G P° and 
p 7^ 0, we have 

e /(*,p) = J- \cp\ 2 e^~^ > coeMti^'fi™ = c e |p|M( ^) > Co e^ M ^ 

P 

for x E K and p 7^ 0, which completes that proof. □ 

Completion of the proof of Theorem li.ifr Recalling that z = e p ^ 2+l9 , we write \^>a N (p)\p 
instead of \<Pa N ( z )\p- By the definition of the Hermitian metric on Lp, we have 



Wi N {p)\ 



\x P g N {p)\p 1 e<*°"> 

"■ vW ""' 'IX^II 2 ' WxU 2 \^Ue p/2 )\ 2N ' 



where, in the right hand side, |$p(e P//2 )| denotes the usual norm in C d+1 . By the definition 
(112)1 of the monomial embedding $p, we have 

\$ p (e p/2 )\ 2 = k(p), P eR m . 

Hence, we have 

-Nf(a N /N,p) p-HfMRniT n) , 

\&M\i = n Yf 112 = ipnp^ . = ^ ww/w - > . 

II Aajy II II Xa N II 

where the function f(x,p) for x 6 P° is defined in (JJOJ). We note that, since a^/N — x = 
0{N- 1 ), we have 

|9 p L ^(x,p)| <C Lj R w (x,p) (48) 

for every multi-index L. By Lemma 13.21 the £ 2 -norm of the un-normalized monomial \a N 
is given by 

||xL|| 2 = ! e- N ^/ N ^detA(p)dp. 



Here it should be noted that detA(p) is a positive integrable function on M m . By Lemma 
13.31 we can choose R > 0, c > such that |p^| < i? and /(ajy/iV, p) > c for any |p| > i? and 
iV. Thus, by choosing a cut-off function g>(p) suitably, we may write 

\\XaJ 2 = e- Nf{x ^ J e- Nb *^R N (x,p)g(p)detA(p)dp + 0(e- cN ). 

Recall that b^(p) = if and only if p = p£, and that p = p^ is the unique critical point of 
6f . The Hessian of 6f at p = p^ is the positive definite symmetric matrix A(P, x). Thus, by 
the Morse lemma, there exists a change of coordinates k from a neighborhood of the origin 
to a neighborhood of p% such that k(0) = pf and that 

£o«(£) = (A(P,z)£,£>/2, |det£)K(0)| = 1. 

By choosing the cut-off function g suitably, we get 

- Nb *^g(p)R N (x,p)detA(p)dp= [ e - N ^ p ^G N (x, £) d£, 



where Gn(x,£) is a compactly supported function in £ such that Gn{x, 0) = det A(P, a;). By 
the derivatives of Gn(x, £) with respect to £ are all bounded uniformly in N. Therefore, 
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by using the Plancherel formula and a formula for the Fourier transform of the Gaussian 
functions, we obtain 

/( N \ ~ m l 2 
e - Nb *^g(p)R N (x,p)detA(p)dp= - V / ditA(^)(l + 0(iV- 1 )), 

which completes the proof. □ 

Remark: Since is a monomial, the asymptotics of is essentially the same cal- 

culation as the asymptotics of the L 2k norm of another monomial. Thus, determining the 
pointwise asymptotics of monomials is equivalent to determining the asymptotics of their 

j2k norms> 

3.2. Pointwise asymptotics: boundary lattice points. Next, we consider the ray Na 
for a lattice point a G P fl Z m , which is allowed to lie in the boundary dP. In such a case, 
we need to work with other coordinates than the usual coordinates on the open orbit (C*) m 
( SZ2J), since the open orbit (C*) m does not cover the set pp 1 (dP). 

In the following, we mean that the faces are disjoint, and the facet is a face of codimension 
one. Thus we call the closed face (or facet) F the closure of F in the minimal affine subspace 
containing F. To describe the coordinates, let t>o be a vertex of P. Since our polytope P is 
Delzant, we can choose lattice points a 1 , . . . , a m in P such that each a 3 is in an edge incident 
to the vertex vq, and the vectors v 3 := a 3 — Vq form a basis of Z m . We choose (open) facets 
Fj, j = 1, . . . , m incident at Vq so that a 3 G" Fj. 

LEMMA 3.4. Let a G P II 7L m , and z G M P . Then, Xai z ) = */ an< ^ on ^V 

p P {z) e{J{F; F is a facet a F}. (49) 

Proof. If z G (C*) m then automatically we have Xa( z ) 0- Thus, we may assume that 
z G /ip 1 (9P). First, assume that pp(z) G F for some closed facet F which does not contain 
a. The last formula in (J23|) for the moment map is globally defined, for the function |m^| 2 
is globally defined. Since pp(z) G F, the coefficients in pp(z) of the lattice points (3 F 
must vanish. Thus, we have Xa( z ) = 0- Conversely, assume that pp(z) G dP is not in the 
set described in (149(1 . In our convention, the faces are disjoint and the boundary BP is the 
disjoint union of faces. Thus, that pp(z) G dP is not in the set described in (j4~9~() is equivalent 
to say that there exists an open face E such that a G E and pp(z) G E. Let Vi,...,Vi be 
the set of vertex of E. Since Vj and a are lattice points, there exists a positive integer uq 
such that n^a = Y^j=i n j v j with rij integer such that J2 n j = n o- Thus we have 

(x:r n °(z) = {x p vl T ni w ®---® (x p vi r ni {z). 

Since pp(z) is in the interior E of the face (polytope) E, each Xvj( z ) can no ^ vanish, and 
hence Xa( z ) 0- ^ 
Using Lemma l3~4l we set 

U vo := {zeM P ; X P o (z)^0}, 
which covers Mp as Vq varies over all vertices. We define 

^(CT-(CT, Ti(z)=ri j (z):=(z*,... t z« n ). (50) 
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The map 77 is a diffeomorphism and the inverse is given by 

z:(CT-(CT, z{r 1 ) = {rf e \...^ m ), 

where e 3 is the standard basis for W 71 (or for C m over C), and T is an m x m-matrix with 
det T = ±1 and integer coefficients defined by 

FV = e 3 , v 3 = a 3 — vq. 

By definition, we have the obvious formula: 

*£(*) = »&(*);£(*), ze(C*) m . 
By Lemma E31 Vj( z ) — > if z G U Vo , z — > fip 1 (Fj). Since a 3 Fj, we have 

The set £/"„„ \ (C*) m is the union of the sets /ip 1 (F J ), and hence the map 77 extends a home- 
omorphism: 

T) : U Vo — > C m , 77(2:0) = 0, z = the fixed point corresponding to Vq. 

By this homeomorphism, the set fXp 1 (Fj) corresponds to the set {77 G C m ; 77^ = 0}. This 
coordinate 77 = (771, . . . , rj m ) is useful to explain toric subvarieties corresponding to faces. 
Namely, let F be a closed face with dimF = m — r which contains Vq. Since Vq G F, we can 
choose Fi 17 . . . , F ir such that F = F^ fl • • • D F r . Then the subvariety /j,p 1 (F) corresponding 
F is expressed, in the coordinate neighborhood U vo , by 

fi P \F) n ^ = {77 G C m ; Tfc. = 0, j = 1, . . . , r}. 

Now, we fix a lattice point a in a (relatively open) face F of dimension dimF = m — r 
such that v G F. Without loss of generality, we may assume that F — Fx fl • • • fl F r . 

To state a result for the lattice point a in the boundary corresponding to Theorem II .21 we 
need to find a function corresponding to the function b^- Since our coordinate 77 is based on 
the lattice points a 3 — vq, it is reasonable to introduce a new polytope defined by the affine 
linear transformation 

f : W m 3 u -> Yu - Tv G R m , 

which maps Z m bijectively onto itself. We set Q := T(P). Then Q is contained in the 
positive orthant {x G lR m ; Xj > 0}, and we have T(Fj) = {x G Q ; Xj = 0}. The face of Q 
corresponding to F is then given by 

Q F :=t(F) = {xeQ;x J =0, j = l,...,r}. 

We denote a point in U Vo = C m as 77 = (£, C) G C m = C x C m ~ r . In this expression, 
( = (0,C) is a local coordinate of the submanifold fXp 1 (F). The modulus square of the 
monomials IxjvcJp with a G P fl Z m is, in this coordinate, given by 



I f(a)|2JV 

~ '' J f-l(T)l 



( 51 ) 

7eQnz m 



20 BERNARD SHIFFMAN, TATSUYA TATE, AND STEVE ZELDITCH 

We then introduce the 'moment map' corresponding to the face F by: 

a u e {p ' u) 
k F (p) 



nm—r 



where a v = |cf_i/ ^ | 2 , and the function k F (p) is given by 

k F (p):= ^e {p ' v) . (53) 

As mentioned above, the submanifold p^}{F) has the coordinate ( i— » (0, £) £ C r x C m ~ r , 
and the torus T m-r acts on it. Thus, it is natural to use the coordinate 

( = e p/2+ie , p,6eW m - r . 

Then, we write 77 = (£, Q = (£, p) for £ = e P//2 . Since we have assumed a e F, we may write 

f(a) = (0,5) G Q F . 

We also define 

Sa(^p) ■= log K(£,p) -(p,a) = \og(k F (p) +£ F (£,p)) - (p, 5), 

Proposition 3.5. In the coordinate r\ = (£,C) a s above, we write \ipx a (£, p)\p for the 
modulus square of the monomial |yjv a |p- We also write T(a) = (0,5), which is in the 
interior of the polytope Q F . Then we have 

where the function \I/ Q (£,p) is given by 

*<*(£, p) = s a {£, p) - s a (0, p F a ) pi = p-/(a) e R m ~ r . (55) 
and (m — r) x (m — r) positive definite matrix A(F, a) is given by 

a e (p£>> 

A(F,a)= > I u®u-a®a. (56) 

We shall prove Proposition 13.51 in the rest of this subsection. First of all, we need the 
following simple lemma: 

Lemma 3.6. In the coordinate (£, £ = e p/2+w ) on C r x (C*) m_r c U Vo! the volume form 
u)™/m\ is given by 

i±i m 1 

-T = —!—^ HZ, P) dm(Odpd9, (57) 

where dm(£) denotes the Lebesgue measure on C r , and the function L(£,p) is given by the 
determinant of the following m x m matrix: 

( a 2 log if a 2 logic \ 

a 2 i% a 2 i% , (58) 

apas, apdp J 
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where K(£,p) is the function defined in ()51jl . 

Proof. For the function /(£,£) on C r x (C*) m ~ r independent of the variable 6 in ( = e p l 2+l6 , 
then the derivatives (df)/(dQ), (df)/(dQ) is given, respectively, by 

d Q Cjdpj' dQ Cjd Pj ' 
We write r\ = (£, £). Then, by this relation, we have 

where L(£,p) is given by (|EHjl. But, we have 

/m-r \ 

<mo = ^ (^n ioi 2 j « 

where dm(Q denotes the Lebesgue measure on C m ~ r . Combining this with (j47j) . we obtain 
the assertion. □ 

The following lemma can be shown by the same argument as in the proof of Lemma 13.11 
and we shall omit the proof. 

Lemma 3.7. The (m — r) x (m — r) matrix defined by 

A F (p) := dp F (p) = v v ® v - p F (p) (8) p F (p) 

is positive definite for every p e IR m ~ r . 

Note that, the map p F : M. m ~ r — > Q° F is a diffeomorphism, and the lattice point a is in 
Q° F . Thus, the vector p^ = /x^ 1 (<5) is well-defined. Hence, the (m — r) x (m — r) matrix 

A(F,a) = A F (pF) 

is positive definite. 

Completion of proof of Proposition HOI The modulus square of the monomial I^tvJp m this 
coordinate is given by 

\XNa(^p)\ P 



IVjVa(f>P)|j 



IIxLII 2 
1 



WxLW 2 

where the function s a (£,p) is defined by (|54jl . Thus, as in the proof of Theorem II. 2\ what 
we need to analyze is the £ 2 -norm 



\XnS = - t [ e- Ns ^L^,p)dm(Odp, (59) 

71 JC r xR m - r 
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where we have used Lemma "3~"6l We note that the function £ F (£,p) defined in (""""f| is of the 
form 

r 

Me,p) = X)A(p)ie*r+r(^p), 

k=1 (60) 

where e£ are the standard basis for R r . The function r(£, p) is of order > 4 in £, and hence 
its derivative up to the second order vanish at £ = 0. Thus, the function s a (£, p) has only 
one critical point (£, p) = (0,p p ). It is not hard to show that the Hessian Hs a (0,p a ) of the 
function ("5""") at the critical point (0, pf ) is given by 

/ 2/i(pg) \ 



Hs a (0,pl 



Vr{ P F a ) 

Mp£) 



2/i (Pp 



V 



(61) 



2/r(pg) 

A(F,a) / 

where the (m — r) x (m — r )-matrix A(F, a) is given by (|56p. and we have used the coordinate 
(xj, t/j, p) with £j = Xj + iyj. Here, it should be noted that the polytope Q = T(P) contains 
the standard basis in M m = W x R m_r . Thus the lattice points (e£, 0) with the standard basis 
e r j in R r is in the polytope Q, and hence the functions fj are all positive. This combined with 
Lemma 13.71 shows that the Hessian Hs a (0,p a ) is positive definite. By the same argument 
as in the proof of Lemma 13*731 the function fcp(p) + @f(£,p) tends to oo as |£| + \p\ — ► oo. 
Therefore, by the standard Laplace method as in the proof of Theorem II .21 we have 

N \-(m+r)/2 e -Ns a (0,pF) 
2^ 



1 

7l r 



y/det Hs a (0,p 



(62) 



L(0,p) 



A direct computation will show that 

f\ det a 
^Hs a ^p a ) = ^ MF 

and hence, the asymptotics (f6~""|) can be written in the form 

1 




det A F (p) 




(63) 



lv p II 2 

lAJVall 



/ at \ -(m+r)/2 

1 — 1 v / det^>)e- 7VsQ(0 ^ ) (l + O(Ar- 1 )). 



□ 



(2vr) r \2tt / 

Dividing |yjva|p by the above, we conclude the assertion 

When our fixed lattice point a is a vertex, say a = v o in the description of the coordinate 
rj, the matrix A(F, a) is not defined suitably. However, clearly the similar asymptotics can 
be deduce by the same method. 
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Proposition 3.8. Suppose that a is a vertex of the polytope P. Then we have 
W P nM\1=( ^Ve-^W-l-l 2 )(l + 0(iV- 1 )) ) 



where the function K(j]) on C m is given by (J5T|) . We also have \c a \ 2 = \x P ( z a)\ 2 where z a 
is the fixed point for the Hamiltonian T m -action such that fip(z a ) = a. 

Proof. In the description of the coordinate rj, we put a = Vq. Then, clearly we have K(0) = 
\ca\ 2 = \Xa( z a)\~ 2 , where the fixed point z a corresponds to the origin 77 = and the function 
K[rj) is defined in floTj) . In this case, we just use the coordinate r\ itself without change of 
variable. In this coordinate, the monomial IxjvcJp ^ s gi ven by 

\x P N M\l = e- m ^\ 

The volume measure Up/m! is of the form: 

up 1 , T , s , , s r , s fd 2 \ogK(ri)\ , . . 

with the Lebesgue measure dm(rj) on C m . In is straight forward to see that the critical point 
of the function log K (rj) is the origin, and the determinant of the Hessian at the origin is 
given by 



H{logK)(0) 



'2L(0) q 
2L(0) 




K(0) 

By using these facts with the Laplace method, we obtain the assertion. □ 

3.3. Moments and C 2k norms. As an application of the pointwise estimates, one can prove 
that eigenfunctions 'localize on tori'. We also determine C 2k norm of the monomials. The 
following is easily shown by using Theorem II. 21 Propositions I3.5l and l3. 81 and the argument is 
the same as in their proofs, and hence we shall omit the proof (see also the proof of Theorem 

EH}. 



Proposition 3.9. (i) Let o>n € NP n Z be a sequence of lattice points satisfying 

for some point x G P°. Then, the measure \<fi p N \pdYo\M P weak* -converges to the 
normalized Haar measure on the m- dimensional torus /ip 1 (a;) ; i.e., 



I vlVaJpdVdMp -> / crd9 

JMp J u^ix) 



for a G C(M P ). 

(ii) Let a be a lattice point in P with dim/ip 1 (o;) = m — r. Then, we have 

W *-j^J^Na\p = d9^- Aa) , 

where dO^-i^ denotes the normalized Haar measure on the (m—r)- dimensional torus 
12 p 1 (a) = T m - r . 

In the above proposition, if a is a vertex, then the left hand side denotes the Dirac measure 
at the fixed point z a G Mp of the Hamiltonian T m -action corresponding to a. 

We note that \(p^ a \ is invariant under the T m action. We denote the Hilbert space of T m 
invariant functions by C 2 nv (M P ). We can restate the conclusion as follows: if a G C°^ V (M P ), 
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then we can regard it as a function on the polytope P. We can also regard |^ Q | 2 as a 
function, say |<p^(/)| 2 , on P, equipped with action variables /. We then have: 

Corollary 3.10. For any a G PnZ m , we have \^ a \ 2 dl — > 5 a ; i.e., J p a\(p^\ 2 dl — > a (a), 
for a G C(P). For G NP D Z m satisfying ()38|) wtt a point x G P°, we /iai>e \<p^ lN \ 2 dI — > 

Next, we determine the asymptotics of the £ 2fc -norm of the £ 2 -normalized monomials. 



Theorem 3.11. (i) Let o>n G NP n Z m be a sequence of lattice points satisfying 

with a point x G P° . Let ||y>~ ||2k denote the L 2k -norm of the C 2 -normalized monomial 
tpP with the weight 7 G NP D Z m . TTien we /jave 

1 / /V \ m/2 1 
H« = ^y (d^j^Tp tl + 0t (iV-)), (64) 

where the Ok{N~ l ) depends on k. 
(ii) Lei a E P be a lattice point with dim/ip 1 (a) = m — r. Then, for r < m — 1, we have 

Il^ = = (^ L J (l + O fc (iV-)). (65) 



a 



For a vertex a, we have 



1 / AT NO- 1 )™ 

IIM = }^(0?J (i + ^iV- 1 ))- (66) 

Proof. The proof of (JMj) and (|66|) is the same as that for (|65|). so we shall give a proof of 
()65|) only. By Theorem 11.21 and Proposition 13.51 we have 



\fcP \\2k 
\YNa\\2k 



kr / AT \ K m + r )/ 2 



(2n) kr ( N 



f e -^*-«.p)L(e, p) dm(0dp(l + OkiN- 1 )), 

Jc r xR m - r 



tt t \2n J (det A(F,a)) k / 2 J C r 

where is given by ([55)1 . As in the proof of Proposition 13.51 the critical point of ty a is only 
the point (0, p^) with = y u^ 1 (a), T(«) = (0,5), and which is non-degenerate. We have 
^0,(0, p^) = 0. Thus, by the standard Laplace method, we have 

||( ~P (2^ /jvV fc - 1)(m+r)/2 1 L(0, Pa ) ^ n t 

H^- 11 *- ^2 J (detA(F,a))^det^ a ( , Pa ) li + Ufcl ' V jj ' 

Therefore, the assertion follows from (J53"j) in the proof of Proposition 13.51 □ 

By using Proposition 13. 5[ we can determine the limit of the sup-norm Hy^jvoJlL- 

Proposition 3.12. (i) Let a N G NP n Z m satisfy flm with a point x G P°. Then we 
have 

, N n m/2 1 

j im h7- ll^ll^o = - 7= = jt^ F - ( 67 ) 

^-00^277/ ™ y / detA(P,x) 
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(ii) Let a be a lattice point with dim /ip 1 (a) — m — r , r < m — 1. Let (£n,Pn) be the 
point where \<PN a \p attains its maximum. Then, we have 

/jvv (m+r)/2 „ p ll2 (27ry , N 

lim — \\<prjl = \> 68 

For a vertex a G P , we have 

lim N- m y P Na \L = \c a \~ 2m = \x*(z a )\ 2 P m , (69) 

where z a G Mp is the unique fixed point for the T m such that fip(z a ) = a. 

Proof. We only give a proof of (J55J) . The function ty a (£,p) defined in (J55)l attains its 
minimum at the point (0, p^), and it tends to oo as |£| + \p\ — > oo. Thus, the function 
(2 ) m I^WeJp i s °f order 0(e~ c7V ) outside a compact neighborhood of (0, p^). On a 
compact neighborhood I? of (0,p£), we have 

lv£a(0,p£)|f.< sup |^ a (C,p)|p< f /A { fL A l + OiN- 1 )). 

K,p)efl V 27r / v/detA(F,a) 

Now ()68|) follows from the above inequality. The same argument with Proposition 13. 51 shows 
PI). □ 

3.4. Asymptotics on projective space. The values of the C 2k norm in Proposition 12.21 
and in the projective-space case of Theorem 11.61 may seem to be different. However, we can 
check that these two coincide by noting the following simple lemma. 

Lemma 3.13. The determinant detA(pS, a) of the matrix A(pS, a) is given by 

(p - H)«i ■■■a m 



det A(pE, a) 



V 



Proof. By applying the differential operator Xjd x . twice to the formula (1 + YaLx x i) p = 
S|/3|<p (p) xl3 ^ or a vec ^ or x £ w e have 

p(p-i)x^(i + ^xo p - 2 +^(i + ^xo p "%= (*Vm-> 

i i \p\< P w 

where 5ij is the Kronecker's delta, and we have set («) = p\/p\(p— We put r = e Mj , 
Then, by (jH3), we have e pQ = and |mjf (e^ /2 )| 2 = e <p - /3) g)/(l + r) p . From this the 
second equation follows. Substituting for x in the above formula, we obtain 



\P\<p W 1 



<P,/3> 1 



We set D m (a%, . . . , a m ) = det A(p, a) with a = (a 1; . . . , a m ). Then a simple computation 
shows that 

D m (a>i, . . . , a m ) £) m _i(a 2 , . . . ,a m ) 1 

= CKj. 

ai • • • a m a 2 • • • a m p 

Thus the lemma follows by induction on the dimension m. □ 
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Note that, for the simplex pE, the faces containing the origin is again a simplex in lower 
dimensional vector space. Therefore, Lemma f3 . 1 31 can be applied to compute det A(F, a). 
In the case of projective monomials, the function 6^, s can be expressed as follows. 

Proposition 3.14. For z e (C*) m , 

b^(z) = plog(l + \z\ 2 ) - log \z a \ 2 + (a, loga) - plogp, 
where log a = (log o?o, . . . , log a m ) . 

This is a direct consequence of the formulas e Pa = ^^a, r = ^e^ Pa ^ and (1 + Yli x i) P = 
J2\/3\<p {p) xlS as mentioned before. 

The pointwise asymptotics of the £ 2 -normalized projective monomials is given in the 
following proposition. 

Proposition 3.15. (1) For z e (C*) m , 

fN\ m/2 v l/2 e -Nb^(z) 

\Zn J a/Go- \a\)ai ■ ■ ■ a m 

uniformly (C*) m . 
(2) Forz = e (p*+u/VN)/2+ie ^ we have 



\&(z)\ 2 



N \rn/2 1/2 -«A(aK«>-i(a,«> 2 )/2 

£) P ' ' (1+0(JV-V)) 



uniformly for \u\ < c, where A(a) is the diagonal matrix with entries a,\, . . . , a m . 

The assertion (1) in the above is a restatement of Theorem 11.21 for the projective mono- 
mials. The assertion (2) follows from (1) and a Taylor expansion of the function b^p. 



4. Asymptotics of distribution functions 



In this section, we find asymptotics of rescaled and un-rescaled distribution functions. Fix 
a lattice point a in P, and let r denote the codimension of the face of P (possibly the open 
face P°) containing a. In analogy with (|4*5jl. we define the constant c(P, a) by 



c(P, a) 



(27T) 



y/det A(P, a) 

(27r) m /<-> i P / \|2\m 



if r < m 

if a is a vertex. 



(70) 



In the following discussion, we give the details for the case where = Na with a lattice 
point a G P fl Z m . For general G iVP fl Z m satisfying (}3*Hj) for a point x G P°, one needs 
only to put r = and replace Na and a by aw and x. 
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4.1. Rescaled distribution functions. We would like to understand the limit distribu- 
tions of the measures \(p^ a \1dVolM P , namely, the limit of its distribution function 

D Na (t) : = Vol Mp {>; | | 2 > *}• (71) 

However, by Theorem II .6[ the k-th moments of the measure \(p^ a \ldVolM P tends to infinity 
as N tends to infinity. Therefore, we need to re-normalize the monomials. 
We write 

/ i\t \ (m+r)/2 / yy\ -(m+r)/4 

dv N=[^) <2Vol Mp , /iva=(^l ¥>£«(*) (72) 

so that 

/ |/U*dv^(z) = l. 

J Mp 

By Propositions l68l and l69| we know that liniAr^oo ||/iVa||L exists and we have 

lim \\f Na \t = c{P,a). (73) 

N—+00 

Furthermore, by Theorem 11.61 we have 

\\r ||2fe /iV\-(m+r)(fc-l)/2 u^p || 2 fc c(P, a) k 1 /Of/vr-l^ f7/A 

ll/ATall^^) = (^j WnJm = k ( m +r)/2 C 1 + °( iV )) ' ( 74 ) 

We consider the limit distribution of the sequence of measures 

"N,r ■= |iWa|*dv^ (75) 
on the real line. The distribution function F^ a (t) of the measure vn, t is given by 

FUt) = {^) {m+r),2 Yo\ Mp {z e M P ; |^L(^)| 2 > {^) {m+r) ' 2 t}. (76) 

The distribution function F^ a defined above can be expressed, in terms of the distribution 
function D^a for the measure |<^J 2 gTVo1m p , as 

T?r u\ (N\( m + r )/ 2 r, ( ( N\( m + r )/ 2 + \ 

We should note that the total mass of the measure vn,t is Vol(Mp)(iV/27r)( m+r )/ 2 , and hence 
it tends to infinity as N goes to infinity. But its k-th. moment satisfies, for each positive 
integer k, 

J = CJ 0]^(1 + O k (N-^)) - C -^$^ (N — > oo). (77) 

By Propositions EH1 and EHl the support of the measure UN, r is contained in a bounded 
interval in [0, +oo) which is independent of N. 

Furthermore, by (|77j) . the measure xdv^^x) is a finite measure on the real line whose 
support lies in a bounded interval in [0, +oo) which is independent of N. This implies that 
the sequence of the finite measures xdi>N tr (x) on the real line has a weak limit, say /i, and it 
must satisfy 

J * k M*) = {k ^Sr»2 > ^ = 0,1,2,.... (78) 
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Thus, the weak limits of the measures xdh>N, r {x) are probability measures, and they have 
supports contained in the interval [0, c(P, a)}. 

Lemma 4.1. Let c be a positive constant, and let h be a positive integer. Let p, be a probability 
measure on the real line such that 

c k 

x k dp(x) = (fc + 1)V2 , k>0. 
Then p must coincide with the measure p c ,h{%) dx where 

PcAx) = ^f^m^XMcAr))^- 1 . (79) 

Proof. We shall use the following formula: 



e t dt 



w s r(s) J 
to compute the Fourier transform 



of the measure p satisfying the condition in the lemma. Substituting s = h/2, w = k in the 
above formula, we get 

k\(k + l) h / 2 T{h/2) 



^ k\(k + l)^ T(h/2)J 



(log(c/ x)f/ 2 - 1 dx, 



cT(h/2) J 

and hence = p2,h(£)- This completes the proof. □ 

As a corollary to Lemma f4. 11 we have the following. 

Corollary 4.2. The sequence of the finite measures {xdu^^^x)} converges weakly to the 
probability measure p^p^^dx, where the density p c (p, a ),r is given by (JTHJ) with c = c(P,a) 
and h = m + r . 

Proof of Theorem \l.J\ The rescaled distribution function F^f a (t) for t > of the measure 
vn = \fNa\ldv r N is given by 



F Na(t) = j X(t,+oo)(x)du Njr (x), 

where X(t,+oo) is the characteristic function of the interval (t, +oo). Now set dpN,r{x) 
xdvN,r(x) and write 

-i 



F Na(t) = J x X(t,+oo)(x)dp N . r (x). 

By Corollary 14.21 we know that the sequence of finite measures p^ converges weakly to 
the probability measure p c (p, a ), r (x) dx. For fixed t > 0, by approximating the function 
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x _1 X( t)+00 )(x) by a sequence of continuous functions, and by using the Lebesgue convergence 
theorem, one easily obtains that, for every < t < c(P, a), 

lim F 7 Na (t)= l x~\u^{x)p^ Pa \ r ix)dx 

N — >+oo / 

c(P,a) 

x- 1 (log(c(P, o^/x))^^ 12 - 1 dx 

log(c(P,a)/t) 

s (m+r)/2-l dg 



c(P,a)T((m + r)/2) 
1 

c(P,a)T((m + r)/2) 
1 



-(log(c(P,a)/t)) (m+r)/2 . 



c(P,a)T((m + r)/2 + 1) 

□ 



4.2. Non-rescaled distribution functions. In the previous section, we derived the limit 
of the rescaled distribution functions F^ a (t). In the definition of the rescaled function F^ a (t), 
the parameter t is rescaled by the factor (iV/27r) (m+r ) //2 so that it tends to +oo as iV — > oo. 
Then the corresponding volume 

Vol Mp (z ; I^LWI 2 > (iV/27r)^/ 2 t) (80) 

is of order N~^ m+r ^ 2 . This was the reason why we need to multiply the volume (}80|) by 
the extra factor (N/2%y m+r ^ 2 in the definition of the rescaled distribution functions. As a 
result, the limit distribution has a universal form (Theorem 11.4)1 . However, one may ask, 
of course, what the limit of the non-rescaled distributions Dj^ a (f) is. But, for each fixed 
t > 0, D Na {t) — > as N — > oo by Theorem II .21 and Proposition 13.51 Therefore, we need to 
replace D^ a {t) by D^ a (WN(t)) for an appropriate sequence {W^it)} of positive functions 
in t > which compensate for the flattening rate. Our sequence will satisfy the conditions 
W N (t) -> as N -> oo and 

Wjv(t) -> 0, Wiv(t)~ 1/Ar -> W(t) (iV -> oo), (81) 

for a function W(t) > 1. We then have the following limit distribution law: 

Theorem 4.3. Let {W N } satisfy JH3J. Then 

lim D N (W N (t)) = Vol Mp ((£,p) e C r x (R*)" 1 -' ; < log W(t)). (82) 

In particular, by taking the function Wn as W^it) = e~ Nt , we have 

lim D N {e~ Nt ) = — [ L{i )P )dm{i)dp. (83) 

//a zs m the interior P° of the polytope P, we have 

lim D N (e- Nt ) = I det A(p) dp. (84) 
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Proof. The proof of flED is the same as that for (JH2J) and (J£3jl . and (JEHJ) follows from (JH2J) 
and Lemma EHO Thus, we give a proof only for (J%2~)) . First of all, we note that the set U Vo 
introduced in Section El is dense in Mp, and hence we may consider the volume of the set 

S N (t) := {(CO e C r x (C*r" r ; |^L(e,C)| 2 > W N (t)}, 

so that £>iv a (Wiv(0) = Vo\ Mp (S N (t)). Let (£, C) e S N . We write C = e p / 2+ie . Then, by 
Proposition 13.51 we have 

*«(£,/>) < log (CiV-( m ^)/ 2 ^(t)(l + a^vl^C)))" 1 ^ , (85) 
where C is a constant, and is a function of order 0(A^ _1 ). Thus, we obtain 
D Na (W N (t)) 

=Vol Mp ((CC) eC'x (C*) m - r ; * a (£,p) < log (CiV-( m+r )/ 2 ^(t)(l + a^(e,C)))^) • 
This combined with Lemma [3.61 shows the assertion. □ 

Our final aim is to prove theorem II .31 which gives the asymptotic limit of the distribution 
function Diy(t) itself without any rescaling. Since, by Theorem 11.21 and Propositions 13. 51 and 
13.81 the monomial |^ a | 2 decays exponentially away from the corresponding invariant torus 
^p 1 (a), it is obvious that -Dtv(^) tends to zero as iV — > oo. So, a problem is to find the decay 
rate of D^it) for any (but fixed) t > 0. 

For every N and < r < m, we set 

s~ = i°s(£) (ro+r)/2 = ^i°s(£). = 

Lemma 4.4. Let t > 0. Let a E P be a lattice point with dim/^p 1 (a) = m — r and lie in a 
face F with dimF = m-r. Let (£ n ,£n) e C r x (C*) m ~ r satisfy |^L(6v,Gv)| 2 > t. We 
write ( N = e O£+ u w)/2+i0iv _ Then, we have 

|6v| 2 + M 2 = 0(r%) = OiN- 1 log iV) 

locally uniformly in t > 0. 

Proof. For every c > 0, we set 

B a (c) = {(£, e ( ^ +u)/2 + i0) eCx (C*) m ~ r ; |£| 2 + \u\ 2 < c 2 }. (86) 

Then, by the argument in the proof of Proposition 13.51 we can find positive constants Cq, Co 
such that 

codei 2 + H 2 ) < * a (£, p F a + u) < Co(|e| 2 + H 2 ) on B a (c). 
Therefore, by Proposition I3.5[ we have 

with some constant C. Therefore, we obtain 

|&v| 2 + M 2 < jjQog(C/t) + s N ), 
which implies the assertion. □ 



DISTRIBUTION LAWS FOR INTEGRABLE EIGENFUNCTIONS 31 

Thus, to obtain the asymptotic estimate of the distribution function D^(t), we need to 
find that of the monomial |^ Q | 2 on the ball B a (crjv) of radius 0{tn) around the point 
(0, Pq), where B a (cr N ) is defined in 



Lemma 4.5. Let c > 0. We denote points B a (cr N ) as (£, £) = (r N w, e^ +TNU ^ 2+ie ) with 
\w\ 2 + \u\ 2 < c 2 . Then we have 

\?Na( r NW,p^ + r N u)\ 2 

=c(P,«)e^-^ < ^ (0 '^ )( ^ ) ' ( ^" ))/2 (l + O(Ar-V2( logiV )3/2 ))) 

where the Hessian H^ a (0, p F ) of the function on C r x M. m ~ r at the point (0, p F ) is given 
by (US} 

Proof. By a Taylor expansion, we have 

y a (r N w, p F a + r N u) = -jl(HV a (0, p F )(w, u), (w, u) ) + R(r N (w, u)) 
for | w | 2 + \u\ 2 < c 2 , where R(r N (w,u)) = O(r^). In particular, we have 

e -N<t> a (r N w,pF+r N u) _ e ~s N ( HV a (0, P Z)(w,u), (w,u) >/2 Q + QfJ^fj.^ 

= e -s N (Hy a (0,pF)( w ,u), (w,u) )/2M _|_ 0(A^ _1/2 (log A^) 3/2 )). 

^From this and the estimate in Proposition IH.5[ the assertion follows for 1 < r < m— 1. For 
r = 0,m, precisely the same argument replacing ^q, by 6^ or logi^ with Theorem 11.21 and 
Proposition IH.8I will show the assertion. □ 

Proof of Theorem M.iA For every t > 0, we set 

MO = to = (£,C) e c x (CT~ r ; |^L(e,C)l 2 > 0- 

Then, by Lemma 14.41 there is a constant c > (depending on a fixed i > 0) such that 
S N (t) C B a (cr N ). Let (£, C) = (rw e^ +rjv ")/2 + iO) G B a (cr N ). Then, by Lemma 1431 
(C) C) £ 'S'iv(i) if and only if 

t < c (P ja yN-s N (H^ a (0,p^( W ,u),( W ,u))/2^ + ajv ( u , jU )) ) 

where a^{w,u) is a function of order A^ _1 / 2 (log iV) 3 / 2 uniformly in (id, m) with |w| 2 +|u| 2 < c 2 , 
and c(P,a) is define in (170)1 . This is equivalent to the following estimate: 

1 fc(Pa) \ 
( HV a {w, u), (w, u) )/2 < 1 + — log -^-^(1 + a N (w, u)) . 

s N \ t J 

Therefore, by Lemma f3. 61 we have 
D Na {t) = —f L(£,p)dm(£)dp, 

71 J {(€,p)=(rNW,p£+r N u) ; ( H(w,u),(w,u) )/2<l+- t - lo S (c(P,a)(l+a(w,u))/t)} 

where we set, for simplicity, H = H^/ a (0, p F ). Changing the variables (£,/?) to we 
have 

r rn+T r 

DNa(t) = - - / L(r N w, p F a + r N u) dm(w)du. 

rt' J{(w,u) ;(H(w,u),(w,u))/2<l+-^- \og(c(P,a)(l+a N (w,u))/t)} 
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A Taylor expansion gives 

L(r N w, p F a + r N u) = L(0, pf) + R N (w, u) 
with the error term Rn(w,u) = 0(tn) uniformly for \w\ 2 + \u\ 2 < c 2 . Thus, we obtain 

r- {m+r) D Na (t) = [ dm(w)du + 0(r N ). 



7P 



{(w,u) ; (H(w,u),(w,u) 



Note that sjy — > oo while r N — > as N — > oo. The function aN(w,u) tends to zero 
uniformly in \w\ 2 + \u\ 2 < c 2 . Therefore, we conclude, for a fixed t > 0, 

lim r N [m+r) D Na (t) = L( ° ,P " ) / dm(u;)du. 



N^oo 7T 



{(w),u)eC r xR m - r ; (H(i«,u),(i«,u)>/2<1} 
m+r 



o(m+r)/2 r/n F\ 

1 ,Pa Wo\ Mp (xeR m+r ; \x\<l) 



7r r v / deT^ 

By f!63|) and the well-known formula for the volume of the unit disk in M m+r , we conclude 
the assertion. □ 
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